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This  paper  presents  an  approach  to  visibility  analysis  for  a  target  located  high  over 
a  terrain  made  of  discrete  elevation  points.  This  analysis  is  a  special  subset  of  a  general 
terrain  intervisibility  analysis  problem.  It  is  an  efficient  algorithm  based  on  a  variation 
of  the  "line-of-sight"  approach,  being  k-times  faster  than  the  ordinary  approach,  where 
k  is  the  ratio  of  the  highest  terrain  elevation  to  the  target  altitude. 


INTRODUCTION 


Given  a  grid  made  up  of  discrete  terrain  elevation  points,  the  visibility  analysis  for 
a  particular  target  point  located  on  or  above  the  terrain  gr  1  is  defined  by  determining 
whether  or  not  a  specific  target  point  can  be  seen  from  all  other  terrain  points. 
Generally,  the  complexity  of  the  visibility  analysis  is  proportional  to  where  N 
denotes  the  number  of  terrain  elevation  points.  'Hiis  complexity  can  also  be  written  as 
a  big-O  notation,  0(N‘  *).  The  complexity  is  derived  from  the  fact  that  target  visibility 
should  be  calculated  for  every  terrain  point  within  the  grid.  This  process  will  take  on 
the  average  the  square  root  of  N  steps  to  calculate.  As  Ae  number  of  terrain  elevation 
points  N  grows  larger,  (i.e.  1.4  million  points  within  an  1®  X  1®  DTED  cell),  it  becomes 
almost  impossible  to  compute  dynamically  or  so  called  "on  the  fly."  For  many  cases, 
however,  this  "on  the  fly"  computation  may  be  desirable  or  required.  For  instance, 
locating  terrain  areas  which  are  hidden  from  a  vehicle  flying  over  the  terrain  is  an 
important  task  of  mission  planning.  Because  these  hidden  areas  would  dynamically  vary 
depending  upon  the  current  location  of  the  flying  object,  they  should  be  calculated  "on 
the  fly"  in  order  to  give  the  mission  planner  the  opportunity  to  try  all  possible  variations 
in  search  of  the  best  route.  It  is  toward  this  "on  the  fly"  conq)utation  requirement  that 
the  scope  of  this  paper  is  focussed  on. 

An  approach  for  the  special  subset  of  the  visibility  analysis  problem  where  the  target 
is  higher  than  the  highest  terrain  elevation  is  given  in  this  paper.  An  efficient  algorithm 
is  proposed  and  proved  as  having  k-times  lower  complexity  than  the  general  algorithm, 
where  k  is  the  ratio  of  the  highest  terrain  elevation  to  the  target  altitude.  In  the 
following  sections,  the  classical  line-of-sight  algorithm  will  be  discussed  first,  followed 
by  a  description  of  the  improved,  faster  version  of  the  algorithm. 
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LINE-OF-SIGHT  ALGORITHM 


The  following  steps  describe  the  algorithm  for  the  visibility  analysis  for  a  target  point 
All  points,  except  for  the  target  point  are  referred  to  as  source  points: 


1.  Visit  an  unvisited  source  point 

2.  Define  a  straight  line-of-sight  line  connecting  the  source  point 
to  the  target  point 

3.  Walk  out  along  this  line  of  sight  from  the  source  towards  the 
target 

4.  If  any  terrain  elevation  point  along  the  line  of  sight  is  higher 
than  the  line  itself,  then  label  the  source  as  "not  visible"  (see 
Figure  1) 

5.  If  the  target  is  reached  without  intersecting  any  terrain  elevation 
points,  label  the  source  as  "visible"  (see  Figure  2). 

6.  Repeat  for  next  source  point 


Figure  1  The  source  point  B  is  not  visible  from  the 
target  point. 


Figure  2  The  source  point  A  is  visible  from  the 
target  point. 


The  above  algorithm  can  be  described  mathematically.  When  a  source  and  a  target  are 
decided  upon,  the  line  of  sight  from  the  source  to  the  target  can  be  expressed  as 
follows: 


X  = 

y  ==  fyit) 
z  ^  fg{t) 


0  S  t  S  1 
0  s  t  s  1 
0  i  t  i  1 


The  terrain  elevation  data  set  is  denoted  by  the  function. 

H  =  G{x,  y) 

Determining  the  visibility  of  the  source  from  the  target  is  equivalent  to  testing  whether 
there  exists  a  t  value  satisfying  the  following  inequality. 

Oitsl 


If  such  a  t  value  exists,  the  source  is  not  visible  from  the  target  The  inequality  test  can 
be  algorithmically  done  by  ranging  the  t  values  from  0  to  1.  The  increment  of  the  t 
value  may  be  determined  by  the  grid  resolution  of  terrain  elevation  points.  Since  in 
most  cases  G  is  a  discrete  function,  interpolation  may  be  inevitable  in  order  to 


approximate  the  left  term  of  the  above  inequality  equation.  The  inequality  test  is  the 
most  time  consuming  process  in  the  line-of-sight  algorithm.  It  is  this  test  process  that 
the  following  proposed  algorithm  seeks  to  improve. 


PROPOSED  ALGORITHM 


The  above  line-of-sight  algorithm  is  still  useful  when  the  target  is  an  object 
located  on  the  terrain.  When  the  target  is  located  above  the  terrain,  the  algorithm  can 
be  significantly  improved  in  terms  of  algorithmic  complexity,  depending  on  the  target 
altitude.  An  improved  vo^ion  of  the  line-of-sight  algorithm  is  proposed  here  for  the 
special  cases  where  a  target  is  located  at  a  much  higher  elevation  than  the  highest 
terrain  elevation  point  (The  highest  terrain  elevation  is  a  priori  in  most  cases.)  Figure 
3  describes  the  intuitive  logic  behind  the  proposed  improvement 


Fifjure  3  Introduction  of  pseudo  target  point 


In  Figure  3,  the  intersection  between  the  line-of-sight  line  and  the  horizontal  plane 
containing  the  highest  elevation  point  is  shown  as  the  pseudo  target  point  The 
possibility  of  any  terrain  elevation  point  crossing  over  the  line-of-sight  occurs  only 
between  the  source  point  and  the  pseudo  target  point  located  on  the  line  of  sight  The 
line  of  sight  is  always  higher  than  the  highest  terrain  elevation  between  the  pseudo 
target  and  the  real  target  For  the  previously  mentioned  mathematical  algorithm,  the  test 
for  the  existence  of  t  values  satis^ing  the  inequality  was  done  by  ranging  t  values  from 
0  to  1.  However,  for  special  cases  where  the  target  is  higher  than  the  highest  terrain 
elevation,  it  is  no  longer  necessary  to  try  all  the  t  values  between  0  and  1.  The  line-of- 
sight  algorithm  needs  to  try  only  the  t  values  for  a  certain  domain,  for  example,  from 


0  to  0. 1  for  cases  where  the  target  altitude  is  10  times  higher  than  the  highest  elevation 
point  The  complexity,  or  the  conputing  time,  of  the  algorithm  is  thus  improved  by  as 
much  as  10  times.  The  walking  out  along  the  line-of-sight  discussed  in  step  5  of  the 
basic  algoritiun  is  no  longer  necessary  beyotxl  the  pseudo  target  point  The  upper  limit 
tp,  of  tile  t  domain  for  which  the  t  values  are  tried  against  the  inequality  test  and  the 
coriesptHiding  pseudo  target  point  (Xp,  yp,  Zp)  am  derived  from  the  previous  line-of-sight 
equations  and  tiie  highest  terrain  elevation,  h,  as  follows: 


f;UJ2) 


=  h 


Accordingly,  Step  5  in  tiie  previous  algorithm  should  now  read  as  follows: 

5.  If  the  pseudo  target  point  is  reached  without  intersecting  any 
terrain  elevation  points,  label  the  source  as  "visible." 


The  previous  inequality  equation  stays  the  same  except  tiiat  the  upper  limit  of  the  t 
domain  now  changes  from  1  to  the  far  lessa  value  of  tp : 

fy(C))  i  fjt),  0  S  t  i  tp 


The  validity  of  this  change  is  obvious,  because  the  terrain  can  not  cross  the  line- 
of-sight  between  the  pseudo  target  and  the  actual  target  as  is  illustrated  in  Figure  3. 
Instead  of  a  formal  proof ,  it  will  be  shown  that  the  complexity  of  the  test  process  is 
upper-bounded  by  a  constant  of  a  value  much  less  than  1.  The  pseudo  target  is  always 
on  the  horizontal  plane  z  s  h  and  varies  only  in  the  x  and  y  coordinates,  depending  on 
the  varying  source  coordinates  (x,^,  y,^,  z,^,^.  Let  (x,^  y^rp*,  z^  be  the  fixed 
target  coordinates.  Then,  the  line  of  sight  will  be: 


4(t) 


“  ^BOOIC9 
~  y eouxc9 

~  ^Bouxem 


^  taxg9t~y eouzcB^ 

*  ^  ^^taxg9t~^ source^ 


The  t  value  corresponding  to  the  pseudo  target,  tp,  is: 


tp  =  = 


h  -  z 


Bource 


^target  ~  ^source 


Since  the  following  inequalities  among  the  highest  terrain  elevation,  the  target  altitude, 
and  the  source  height  always  hold 


^Bouxee 


tBxgat 


^  h  <  z, 


the  tp  is  upper-bounded  as  follows: 


t  =  ^  ~  ^aourc0  ^  h 

^  ^targmt  ~  ^bouzc*  ^tazget 

The  highest  elevation  h  and  the  target  altitude  z,,^  are  known  constant  values. 
Therefore,  tp  is  upper-bounded  by  a  constant 


This  makes  the  complexity  of  the  test  upper-bounded,  and  the  proposed  algorithm 
becomes  k-times  less  complex  than  die  classical  algorithm  with  t  domains  always 
between  0  and  1  (where  k  is  the  reciprocal  of  the  upper-bound,  which  is  the  ratio  of  the 
highest  terrain  elevation  to  the  target  altitude). 


CONCLUSIONS 

This  paper  has  shown  that  the  classical  line-of-sight  algorithm  can  be 
significantly  improved  for  special  cases  where  a  target  is  higher  than  the  highest  terrain 
elevation.  The  algorithmic  complexity  of  the  improved  version  has  been  analyzed  and 
shown  as  having  an  upper-bound  which  is  quanted  by  the  ratio  of  the  highest  terrain 
elevation  to  the  target  altitude. 


